On the level density of spin chains of Haldane— Shastry type 
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We provide a rigorous proof of the fact that the level density of all su(m) spin chains of Haldane- 
Shastry type associated with the An-i root system approaches a Gaussian distribution as the 
number of spins A'' tends to infinity. Our approach is based on the study of the large TV limit of the 
characteristic function of the level density, using the description of the spectrum in terms of motifs 
and the asymptotic behavior of the dispersion relation. 
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A well-known conjecture on spin chains of Haldane- 
Shastry (HS) type [MW^ states that their level density 
becomes Gaussian as the number of sites tends to infin- 
ity. Although this conjecture has been numerically ver- 
ified for all chains of HS type whose spectrum has been 
computed in closed form [^, M, lllHlJl , a rigorous proof 
thereof is lacking, except in a few exceptional cases in 
which the partition function factorizes [l8|. In this pa- 
per we settle the conjecture in the aflfirmative for spin 
chains of HS type associated with the An-i root system, 
which is both the simplest and the most studied case. 

Our result has implications in connection with two fun- 
damental conjectures in the theory of quantum chaos 
that we shall now discuss. The first of these conjectures, 
due to Berry and Tabor [1^, asserts that the proba- 
bility density of spacings between consecutive levels in 
the spectrum of a quantum system whose classical ana- 
log is integrable follows Poisson's law p{s) = e*. The 
second conjecture, formulated by Bohigas, Giannoni and 
Schmidt [20| , posits that for a fully chaotic quantum sys- 
tem this density is instead given by Wigner's surmise 
p{s) = (7rs/2) exp(— 7rs^/4), characteristic of the Gaus- 
sian orthogonal ensemble (GOE) in random matrix the- 
ory [2l|. It is important to bear in mind that in this 
context the term "spectrum" refers to what is known as 
the unfolded spectrum, which by construction has an ap- 
proximately uniform level density. The energies of this 
unfolded spectrum are obtained from the "raw" energies 
Ei through the mapping Ei i-^ r]i = rj{Ei), where f] is 
a continuous approximation to the cumulative level den- 
sity. Thus, the knowledge of this continuous approxima- 
tion is essential for testing the latter conjectures. It turns 
out that in all spin chains of HS type studied so far, if 
one assumes that the continuous part of the cumulative 
level density is Gaussian, i.e.. 
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of the spectrum) , the spacings density follows a charac- 
teristic distribution which is neither of Poisson's nor of 
Wigner's type [§,[ill-[l3- 

Let us briefiy recall the definition of the three spin 
chains we shall deal with in what follows. The Hamil- 
tonian of the original sn(m) Haldane-Shastry chain is 
defined as 
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where (as always hereafter, unless otherwise stated) 
the sum runs from 1 to the number of spins N, and 
e — 1 (resp. e — —1) for the ferromagnetic (resp. 
antiferromagnetic) chain. The operator Sij permutes 
the ith and jth spins, i.e., its action on an ele- 
ment |si, . . . ,Si,.. .,Sj,.. .,sn) (sfc e {1, . . . ,m}) of the 
canonical spin basis is given by 
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The permutation operators Sij can be expressed in terms 
of the generators f^ of the fundamental representation of 
su(m) at the fc-th site as 



s., = 2j2 nt'; 
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with the normahzation tr(i^t|) = \5°-^ . The chain ((H) 
is intimately connected with the Hubbard model. For 
instance, it can be obtained from the one-dimensional 
Hubbard model with long range hopping introduced in 
Ref. [24I when the on-site interaction tends to infinity 
and the sites are half-filled. The rational version of 
the HS chain ([T]) was subsequently introduced by Poly- 
chronakos 3] and Frahm 23]. The Hamiltonian of the 
Polychronakos-Frahm (PF) chain can be taken as 



H 



y- 



eS,. 



i<j 



i^^-Q' 



(2) 



•Electronic address: 'aenciso@fis.uc m.esl 

tElectronic address: ffinkel@fis.ucm.es] 

t Corresponding author. Electronic address: lartemio@fis.ucm.esl 



where the chain sites ^fe are no longer equidistant, but are 
given by the zeros of the Hermite polynomial of degree N. 
Finally, the hyperbolic version of the HS chain, known as 



the Frahm-Inozemstsev (FI) chain [4!], is defined by the 
Haniiltonian 
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The chain sites in this case are given by ^i = ^ log d , 
where Q is the ith zero of the Laguerre polynomial L'^ 
with a > 0. In particular, unlike the previous two chains, 
the sites of the FI chain depend on an essential parame- 
ter. 

Each of the chains (H])-® can be obtained from 
a corresponding spin dynamical model of Calogero- 
Sutherland type p4i-26i] by applying the so-called freez- 
ing trick [3]- As first shown by Polychronakos [23|, this 
connection can be exploited to derive closed-form expres- 
sions for the partition functions of the above chains. It 
turns out that these expressions can be rewritten in a 
remarkable unified way as [HI [11, [IE HI 

r N—r 

where q = g^i/C^'B^)^ k ^ (j^^^ _ ^ f^^^ jg an element of the 
set Vn of partitions of N with order taken into account, 
and the spin degeneracy factor d{ki) is given by 
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The numbers Ki in Eq. ^ are defined as Ki — 
E]=i kj e {1, . . . , A^ - 1}, and {K[, . . . , K'^_^} = 
{1, . . . , A^ — 1} \ {Ki^ . . . , Kr-i}- The partition function 
Z depends on the chain under consideration only through 
its dispersion relation J-"(i), given by 

i{N — i) , for the HS chain 

T{i) = {i, for the PF chain (5) 

i{a + i — \) , for the FI chain . 

Using Eq. gD, Basu-MaUick et al. [H HI derived a 
simple set of rules for generating the spectrum of the 
chains ([I])-© in terms of Young tableaux of certain ir- 
reducible representations of the Yangian y(gl(m)). It 
can be easily shown that these rules are equivalent to the 
explicit formula 

N-l 

En^^S{ni,n.i+i)T{i) , n = (ni, . . . , tiat) , (6) 

where the quantum numbers Ui independently take the 
values 1, . . . , m. As to the function S, it is given by 



Hj, k) 



1 , i <k 
, J ^ fc , 



(7) 



in the ferromagnetic case, whereas in the antiferromag- 
netic one it suffices to exchange and 1 in Eq. ^. 
The vectors ^(n) S {0, 1}^~^ with components (5i(n) = 
(5(ni, Tii-i^i) are in fact the celebrated motifs introduced 
by Haldane et al. in Rcf. [30,]. It should be emphasized 
that the formula (O for the energies is obtained from the 
partition function and not vice versa, as is usually the 
case. 

Equation ^ shall be our starting point for establishing 
the asymptotically Gaussian character of the level den- 
sity of the chains ([I])-®. In fact, since the sum of the 
Hamiltonians of the ferromagnetic and antiferromagnetic 
chains is a constant, from now on we shall restrict our- 
selves to the ferromagnetic case unless otherwise stated. 
We shall start by deriving a unified formula for the mean 
/i and variance a^ of the spectrum of the latter chains 
in terms of the dispersion function F. Consider first the 
mean energy, defined by 
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By Eq. ([7]), the coefficient of J"(z) in the previous expres- 
sion equals 1 for Ui+i > Ui, regardless of the values taken 
by the N — 2 remaining quantum numbers n^., and is 
otherwise. Thus 
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Similarly, the variance of the energy is given by a^ 
(E^) — ^^, where 
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Taking into account that df — Si and proceeding as before 
we easily obtain 
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Hi < j — 1, the coefficient of F{i)F{j) in the last sum 
equals 1 provided that Ui < rii+i and Uj < Uj^i, and is 
otherwise zero. Likewise, the coefficient of J^{j — l)T{j) 
(j = 2, . . . , A'^— 1) is 1 if JT.j-1 < Uj < nj+i, and vanishes 



otherwise. Hence 
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Using Eq. ((§]) for ^, after some straightforward algebra 
we obtain 
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Since, up to an additive constant, the energies of the 
antiferromagnetic chains (II|)-(131) differ from those of their 
ferromagnetic counterparts by a sign change, it is clear 
that Eq. ^ is also valid in the antiferromagnetic case. As 
to Eq. ([5]), using the antiferromagnetic analog of Eq. ([7]) 
and reasoning as before it is immediate to show that the 
mean energy of the antiferromagnetic chains is given by 
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It may be easily verified that the unified expressions ^- 
(IIOP coincide with the values of /i and a^ previously com- 
puted on a cas e by case basis for the (antiferromagnetic) 
chains (HD-® [ll[li,[l3. 

After these preliminaries, we are now ready to present 
the main part of our proof. As in our previous paper [lg| , 
the proof is based on analyzing the limit as N tends to 
infinity of the (normalized) characteristic function of the 
level density, defined by 

m 

rii ,...,niv — 1 

= m-^e-*^*/'^Z(e**/'^). (11) 

Note that 0{t) is simply the Fourier transform of the level 
density, after normalizing the spectrum to zero mean and 
unit variance. The importance of the characteristic func- 
tion in the present context lies in the following standard 
result (see, e.g., Ref. [SI]): in the limit N -^ oo, the 
level density (normalized to unity) approaches a Gaus- 
sian with parameters fi and a provided that 
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for all real t. In order to analyze the asymptotic behavior 
of the characteristic function pT|) . we first rewrite the 



partition function Z{q) of the chains (HJ-® using the 
explicit formula ^ for the energies, obtaining 
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where the transfer matrix Tj{q) (j — 1, . . . , N — 1) is the 
m X m matrix with elements 



1 
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Note that, since 6{k, I) depends only on the difference k—l 
(cf. Eq. (m), Tj{q) is a Toeplitz matrix. More precisely, 
using the explicit definition ^ we see that 



where the matrix T{uj) is given by 
/I w™ ••• a;™\ 
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and 
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The m X m matrix (|15p can be easily diagonalized for 
arbitrary m and to. Indeed, consider the vectors v''{(jj) 
(k = 1, . . . , m) with components 

.,^(c.)=[c.e^-'=/™]"^^ j = l,...,m; 



note, in particular, that v"^{uj) = {ui^' 
then have 



m[r(c.)t;'=(c.)], 
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where 
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Thus v^{uj) is an eigenvector of T{oj) with eigenvalue 
Aa:(w). The above result is valid for arbitrary w € C. 



When u is unimodular, the vectors v^{oj)/^/rn {k = 
1, . . . , m) form an orthonormal basis of C". Indeed, they 
are clearly of unit length, and their scalar product is given 
by 



v'^-v'^ 
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where we have used the fact that uJ = llj~^. In other 
words, when Iwl = 1 we can write 



T{iu) = U{Lj)D{uj)U''{Lj), 



(18) 



where 



I?H = diag(AiH,...,A„,H) (19) 

and U{uj) is the unitary m x m matrix with entries 
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Let us now go back to the characteristic function pTjl . 
Since 



uj,(^'/^)^^^\ 7, 



_ Hj) 



(22) 



we can apply Eqs. dUll-dlDl) to the matrices Tjie'*/'^) = 
r(e*'^'*) in Eq. 1^. We thus readily obtain 



ip{t) = V M„„'(i) 



(23) 



n,n' — 1 

where the m x m matrix M{t) is given by 

M(i) = t/(e*^i*)i:»(e*''i*)Bi(i)---i:»(e*''"-^*)Sjv-2(i) 
X i:)(e'T«-i*)[/t(e*''"-i*) (24) 

with 

B^{t) = C/t(e''^^*)C/(e'''^+i*). (25) 



From Eqs. ([23)l - (l25l) one can determine the large N limit 
of the characteristic function (p{t), as we shall next dis- 
cuss. To this end, we note first of all that 



7, = OiN-'/') 



(26) 



Indeed, for the HS and FI chains J"(j) is 0{N'^) by 
Eq. (0), while it can be checked that a is 0{N^/^) by sub- 
stituting the corresponding expressions of J-' in Eq. (0) 
into Eq. ^ (cf. the explicit formulas in Refs. 11] and 
[l3|). Likewise, for the PF chain J^{j) is 0(iV), while a 
is 0(Ar3/2) [i||. By Eq. (gOD, this implies that 



where 



U{e'^'') ^R + 0{N-^/^). 



R=U{1) 
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(28) 



is a constant unitary matrix (independent of N) . In order 
to estimate Bj(t), note first that for all three chains H])- 
dSl) we have 



7j-7j+i-0(iV-3/2); 



(29) 



indeed, J"(j) - J'(j + 1) is 0{N) (rcsp. 0(1)) for the HS 
and FI (resp. PF) chains. Taking this into account and 
using Eq. (|20p we immediately obtain 
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so that 

Bj(t) = I + 0(7V-3/2). (30) 

On the other hand, from Eq. ([TT]) it easily follows that 

|Afc(e'^^*)Kl, A: = l,...,m. (31) 

Equation dM} and the estimates (H?!), ([201) and dSII) im- 
mediately yield the asymptotic formula 

M(t) = i?A(t)i?^ + 0(iV-i/2)^ (32) 

where A(t) is the diagonal matrix with entries 



Afe(i)= n^fc(e'^^*)' k = \,..., 

3 = i 



(33) 



Inserting Eqs. ([5^ -([55 1) into Eq. (P5| and using the iden- 
tity ()2ip with w = 1 we obtain the simple asymptotic 
estimate 



0{t) = e-'''*/"A™(t) + 0{N-^''') . 



(34) 



In view of the latter equation, in order to complete our 
proof of Eq. (|12p we just have to determine the asymp- 
totic behavior a.s N ^ oo oi the eigenvalue Am(w), with 
ui = e*''^* unimodular. By Eq. pT|) . this eigenvalue is 
given by 
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provided that t ^ {2tt/^j) Z. Note that for any fixed t ^ 
the condition mt ^ {Tj-k/^j ) Z is fulfiled for sufficiently 
large N on account of Eq. (P5)l . and implies the weaker 
condition t ^ (27r/7j)Z. Thus for all t ^ Eq. f5^ 
holds if iV is large enough, and Xm{e-'^^^^) 7^ 0. The latter 
equation, together with the elementary Taylor expansion 



log 
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mie^^ — 1) 
and the identity 



= \{m-l)ix-^^{m'-l)x^ + 0{x') 



easily yields the asymptotic formula 
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In order to estimate the coefficient of t^ in the previous 
formula, it suffices to note that Eq. ^ can be equiva- 
lently written as 



Af-l 
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where we have used Eqs. (P5| and 
mate. Hence 



for the last esti- 



_!^+logA„(t) = -^ + 0(A^-i), 

(7 2 

which obviously implies Eq. (|12p in view of Eq. (|34p . 

We shall conclude by summarizing the main result of 
this paper and presenting an outline of related future 



work. We have rigorously shown that for all spin chains 
of Haldane-Shastry type associated with the A^-i root 
system the level density (normalized to unity) approaches 
a Gaussian distribution as the number of sites tends to 
infinity. Our proof essentially relies on two key properties 
of these chains, namely Eq. ^ for the energies in terms of 
the motifs ^, and the estimates (|^5|) and (|29p involving 
the large N behavior of the dispersion relation ^ . 

Our results admit several natural generalizations. For 
instance, one could consider the su(n|n') supersymmetric 
extensions of the chains (IT])-(l3]), some of which have al- 
ready been studied in the literature [i3,[i3l ■ It is straight- 
forward to check that Eq. ([7]) for the motifs should be 
replaced by 



5{j,k) 



j>k or j = k>n 
j < k or j — k ^ n , 



where j,k = 1, . . . ,n + n' = m. As a consequence, the 
last n' elements in the main diagonal of the transfer ma- 
trix (fTS]) are replaced by w™ , so that the resulting matrix 
is no longer Toeplitz. Although this fact certainly compli- 
cates the explicit diagonalization of the transfer matrix, 
we believe that the main ideas behind our proof can still 
be applied to this case. 

It is also natural to consider the generalization of our 
result to spin chains of HS type associated with other 
root systems, like BCn or Dj^. The main difficulty in 
this respect is the fact that for these chains no description 
of the energies in terms of motifs akin to Eq. ^ is known 
so far. At least for the Sutherland spin chain of BCn 
type [5|,|8|, some preliminary results of our group indicate 
that such a description is possible, and that our proof can 
be suitably adapted to this case. 
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